#include<iostream>
#include"io.h"
#include"conjugate_gradient.h"

using namespace math;

int main()
{
	matrix<double> A(3, 3, 1);
	for(int i = 0; i < 3; ++i)
		A(i, i) += 3;
	output(A);
	std::cout << "\n";

	vector<double> b(3,1);
	output(b);
	std::cout << "\n";

	auto x = cg_solver(A, b, b);
	output(x);
	std::cout << "\n";
	output(prod(A, x));
}
